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Abstract 

Exclusive semileptonic decays B Z?i(2420)£^' and B D2{2460)ii' are 
studied at the subleading order of the heavy quark expansion. The sub- 
leading Isgur-Wise functions resulted from the kinetic energy and chromo- 
magnetic corrections to the HQET Lagrangian are calculated by QCD sum 
rules in the framework of the heavy quark effective theory. The decay rates 
and branching ratios are computed with the inclusion of the order of 1/ mq 
corrections. It is found that the l/mg correction to the decay rate is not 
large for B D2 but is very large for B ^ Di. 

PACS number(s): 14.40.-n, 11.55.Hx, 12.38.Lg, 12.39.Hg 



Typeset using REVTgX 



I. INTRODUCTION 



In recent years there has been a continuous interest in the investigation of semileptonic 
decays of B meson into excited charmed mesons. This interest arises from several reasons. 
The current experimental data show that the exclusive B transitions to the ground state 
s-wave D and D* mesons make up only approximately 60% of the inclusive semileptonic 
decay rate, thus a sizeable part of semileptonic B decays should go to excited D meson 
states. Indeed, these decays have been observed and more experimental data are collected 
with an increasing accuracy Theoretically, the semileptonic B decays into excited 

charmed meson states can provide an additional source of information for determining the 
CKM matrix element Vcb as well as exploring the internal dynamics of systems containing 
heavy-light quarks. 

The heavy quark symmetry has important consequences on the spectroscopy and 
weak decay matrix elements of mesons containing a single heavy quark Q. In the infinite 
mass limit, the spin and parity of the heavy quark and that of the light degrees of free- 
dom are separately conserved. This allows that the hadronic states can be classified in 
degenerate doublets by the total angular momentum j and the angular momentum of the 
light degrees of freedom j^. In the case of qQ mesons, coupling with the spin of heavy 
quark sq = 1/2 yields a doublet with total spin j = je ± 1/2. The ground state mesons 
with j[ = I are the doublet {D,D*) for Q = c and {B,B*) for Q = b. The excited 
heavy mesons with jf = 1/2+ and 3/2+ can be classified in two doublets of spin symmetry 
(0+,l+) and (l+,2+), which are identified as [D^, D[) and {Di, D2) for charmed mesons, re- 
spectively. The other important application of heavy quark symmetries has been the study 
of semileptonic transitions between two heavy hadrons. The hadronic matrix elements of 
weak currents between members of the doublets identified by ji and ji' can be expressed 
in terms of universal form factors which are functions of the dot-product, y = v ■ v', of 
the initial and final hadron four-velocities. A well-known result is that the semileptonic 
B decays to ground state D^*^ mesons, in the mg 00 limit, can be described in terms 
of a single universal function, the Isgur-Wise function ^{y). While the decays to P-wave 
excited D mesons with j[ = 1/2+ and 3/2+ require two independent functions, ({y) and 
T(y) [0,^, respectively, in the limit rnq — > 00. 

There are Aqcd/'^q corrections to the weak matrix elements parameterized by form 
factors at the ttiq — » 00 limit. The Aqcd/'h^q corrections to the leading term can be 
analyzed in a systematic way in heavy quark effective theory (HQET) in terms of a reduced 
number of universal parameters. The Aqcd/'^q corrections may play an important role 
for i?-decay modes into excited charmed states since the corresponding transition matrix 
elements in the infinite mass limit vanish at zero recoil point because of heavy quark spin 
symmetry, while Kqcd/^q corrections to these decay matrix elements can give nonzero 
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contributions at zero recoil ^j. The kinematically allowed range for these decays mostly 
occurs near the zero recoil point, thus the magnitude of Kqcnl'mQ corrections might be 
comparable with the leading order result. 

The universal functions must be estimated in some nonperturbative approaches. A 
viable approach is the QCD sum rules formulated in the framework of HQET |jlO|. This 
method allows us to relate hadronic observables to QCD parameters via the operator prod- 
uct expansion (OPE) of the correlator. A fruitful application of QCD sum rules has been 
the determination of the Isgur-Wise functions parameterizing the B ^ D^*^ semileptonic 
transitions up to the Aqcd/'^q corrections |]TT]-^]. The QCD sum rule analysis for the 
semileptonic B decays to excited D mesons involves the determination of the universal 
form factors. At leading order in the l/friQ expansion, the two independent universal form 
factors ({y) and T{y), that parameterize transitions B — >■ D** [D** being the generic L = 1 
charmed state), have been calculated with QCD sum rules ||15| , p!6[| . Moreover, perturbative 
corrections to 0{as) have been included in the QCD sum rule for C{y) in |T^. The other 
approaches include various versions of the constituent quark model |]TB|-Q and relativistic 
Bethe-Salpeter equations [^]. The analysis of Aqcd/^^q corrections is an important issue 
for the semileptonic B decays to excited D mesons. Such corrections have been investigated 
in terms of meson mass splittings in Ref. [H and by employing the relativistic quark model 



in Ref. |26]. The corrections have also been included by a variant approach in HQEFT 



in Ref. At the order l/mg, the corrections for matrix elements oi B ^ D** include 
contributions from higher- dimensional operators in the effective currents and in the effec- 
tive Lagrangian. For the semileptonic transitions B Diiu and B -D2^z^, the former 
give rise to two independent universal functions, denoted by Ti{y) and T2{y) In the 
framework of QCD sum rules, these two independent form factors have been investigated 
in our previous work in Here we shall focus on the second type of corrections, which 
originate from higher-order HQET effective Lagrangian. 

The remainder of this paper is organized as follows. In Sec. y we review the formulas 
for the matrix elements of the weak currents including the structure of the Aqcd / "^q cor- 
rections in the heavy quark effective theory. The QCD sum rule analysis for the subleading 
Isgur-Wise functions related to the corrections from the insertions of the kinetic energy 



and chromomagnetic operators is presented in Sec. |ITI| . Section |^ is devoted to numerical 
results. Concluding remarks are given in Sec. 0. 



II. THE HEAVY-QUARK EXPANSION AND THE SUBLEADING 
ISGUR-WISE FORM FACTORS 

The theoretical description of semileptonic decays involves the matrix elements of vector 
and axial vector currents [V^ = cj'^b and A^^ = cj^j^b) between B mesons and excited 
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D mesons. For the processes B 
parameterized as 



Diiv and B D^f-i', these matrix elements can be 



(DiK, e)| \B{v)) = Me*^ + {fy,v^ + fv,v'^) e* ■ v , 
{D,{v\ e) I \B{v)) = I U e^'^^^elvpv'^ , 
{D;iv', e)\A'^ \B{v)) = kA, e*^''v^ + {kA^ + ^a^v'^ e^p v^^v" 
{Dl{v',e)\ \B{v)) = t kye^'^P'^eiyvpv'^ . 



(la) 
(lb) 
(Ic) 
(Id) 



Here form factors and ki are dimensionless functions of y. In the above equations we 
have used the mass-independent normahzation {M{v')\M{v)) = (2'kY2p^ /rriMS^ {p — p') for 
the heavy meson states of momentum p = ttlmv. Therefore, there is a different factor from 
the corresponding equations in Ref. [8]. The differential decay rates expressed in terms of 
the form factors are given by (taking the mass of the final lepton to zero) 



dr^, _ Gl\V,b\^m% 
dy 



+ 



487r3 



{y-ri)fv, + {y'-i){fv, + nfv,) 



dTn* _ Gl \Vcb\'^m% 3 



dy 144773 



l\2{l-2yr, + rl) f^^ + (y'-l)fl 



{y' - 1) kl 



ri (?/2- 1)3/2 3(1- 2|/r2 + 



+2 {y - ra) /ca^ + {y^ - 1) (^^3 + ^2 kA^ 



(2) 



(3) 



where ri = mD-^/rriB and r2 = m£)*/mB- 

The form factors fi and ki can be expressed by a set of Isgur-Wise functions at each 
order in AqcD/'^c.fe- This is achieved by evaluating the matrix elements of the effective 
current operators arising from the HQET expansion of the weak currents. A convenient 
way to evaluate hadronic matrix elements is by using the trace formalism developed in 
Ref. p9| to parameterize the matrix elements in Eqs. ([^). Following this method, one 
introduces the matrix representations 

1+1^ 



Hv 



2 

1 + 



37.(7^ 



(4a) 
(4b) 



where P^, P*^ and , P*^^ are annihilation operators for members of the jf = 1/2 and 
3/2+ doublets with four-velocity v in HQET. The matrices H and F satisfy '^H^ = = 
-H,^, i,Fi^ = Fi^ = -Fi^l FI^^^ = 0,andVf,F,^ = 0. 

At the leading order of the heavy quark expansion the hadronic matrix elements of 
weak current between the states annihilated by the fields in and F^, are written as 

h^^^ T h^^ = r Tr [v^F^, V H,] . (5) 
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where h^^^ is the heavy quark field in the effective theory and r is a universal Isgur-Wise 
function of y. 

At the order Aqcd / ""^q there are contributions to the decay matrix elements originating 
from corrections to the HQET Lagrangian of the same order 



5C 



2mQ 



(6) 



where 



The matrix elements of Aqcd/'^q corrections from the insertions of the kinetic energy 
operator Okin and chromomagnetic operator Omag can be parameterized as 



d^xT{0S,J 



X) 



X] 



d^x 



' 1 ^mag,D ' 



X] 



/^i'^r/^W](0)} = r/£?Tr{^.F,^r/7.}, 
V h^^] (0) } = Tr [v^F^, V H^] ; 

;/^?r/^l^)](o)} 
;/^i'^r/.(^)](o)} 



Tr { n% F:, 



0/3 



1+/ 



(7) 
(8) 

(9) 
(10) 



The functions i]^Q^\y) have mass dimension and effectively correct the leading order Isgur- 
Wise function r(?/) since the kinetic energy operator does not violate heavy quark spin 
symmetry. The most general decomposition for TZ'^^'^'> are 



'n(c) (c) I (c) I (c) 

-oW (^) , W ' I (^) ' 



(11) 



where rji are function of and have mass dimension one. 

There are also order Aqcd / ^c,b corrections originating from the matching of the b ^ c 
flavor changing current onto the effective theory, they can be parameterized in terms of 
two independent Isgur-Wise functions, ti and T2 

Summing up all the contributions up to order Aqcd/^^c,;;, it is straightforward to express 
the form factors fi and ki parameterizing B ^ Diiv and i? — D2 £ z/g semileptonic decays 
in terms of Isgur-Wise functions. The explicit expressions for fi and ki are as follows p|: 

V6 /a = -{y + l)r - eb{{y - 1)[(A' + A)r - {2y + l)ri - r^] + + l)r/,} 
-£,[4(yA' - A)r - 3(y - l)(ri - r^) + {y + l)(%e - Sr^i - 87/3)] , 

v^/y, = (1 - y^)T - e,{y^ - 1)[(A' + A)r - (2|/ + l)n - + r^^] 

-£,[4(y + l)(yA' - A)r - _ l)(3ri - 8x2 - r^ke + 27^1 + 81/3)] , 
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V6fv, = -St - 3£fe[(A' + A)r - {2y + l)n - + Vb] 

- l)ri + 5t2 + 3r/ke + lOr^i + A{y - l)r/2 - Sr^g] , 
V6fv, = {y- 2)r + eb{{2 + y)[{A' + A)r - {2y + l)n - r^] - (2 - y)r],} 
+e,[4{yA' - A)r + (2 + y)Ti + (2 + 3y)r2 + - 2)r/ke 

-2(6 + 2/)r7i - 4(y - l)r/2 - (3y - 2)7/3] (12) 
= - + - (2y + l)n -T2 + r]i,] - Ecin -T2 + ?7kc " 2?7i + 773), 

kA, =-{l+y)r- et{{y - 1)[(A' + A)r - (2y + l)n - r^] + (1 + i/)r/,} 
-e,[{y - l)(ri - rs) + {y + l)(r/ke - 2r7i + r]^)], 
= -^edn + r]2), 

kA3=T + EbiiA' + A)r - {2y + l)ri - ra + r]b\ - edn + T2 - ?7kc + 2?7i - 2?72 - 773) , 

where Eq = l/(2mQ), r/ke = rjL and rjb = r]^J + Qiji^ - 2{y - 1) 7]^^'' + r]i^\ A(A') is mass 
parameter of ground state (excited) mesons in HQET and the superscript on r- ^'^ and rjl'^^ 
are dropped. 

The form factors r and Tj {i = 1, 2) in HQET, that occur in Eq. (|T^) have been in- 
vestigated by using QCD sum rules in our previous work ||T6| , |28[1 . In the following sections 
we shall extend the QCD sum rule analysis to the calculation of the subleading Isgur-Wise 
functions, f]ke{y) and rii{y), associated with the insertions of kinetic energy and chromo- 
magnetic operators of the HQET Lagrangian, SC in Eq. (0). 



III. QCD SUM RULES FOR ISGUR-WISE FUNCTIONS rjKE AND yi 

A basic element in the application of QCD sum rules to problems involving excited 
heavy mesons is to choose a set of appropriate interpolating currents in terms of quark 
fields each of which creates (annihilates) an excited state of the heavy meson with definite 
quantum numbers j, P, je. The proper interpolating current Jj'^pj"^ for the state with 
arbitrary quantum number j, P, ji in HQET was given in These currents have nice 
properties. They were proven to satisfy the following conditions 

(OKKr(O)b-'' P'^J^) = fpJn'Spp'S^a'f'"'"' ' (13) 
z (0|T (jS? (a:)jjfp;f(0)) |0) = 5,,v5pp,5,,,;(-l)^' 5 ^^^^ ■ ■ ■ g^^'^ 

X J dt5{x-vt)Upjdx) (14) 

in the ttlq — >■ 00 limit, where 7]°'^"'°'^ is the polarization tensor for the spin j state, v is the 
velocity of the heavy quark, g"^^ = g""^ — v°'v^ is the transverse metric tensor, S denotes 
symmetrizing the indices and subtracting the trace terms separately in the sets (ai ■ ■ - aj) 
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and [Pi - ■ ■ Pj), fpj^ and Hpj^ are a constant and a function of x respectively which depend 
only on P and je. Because of Eqs. (|13D and (0), the sum rules in HQET for decay 
amplitudes derived from a correlator containing such currents receive a contribution only 
from one of the two states with the same spin-parity (j, P) but different je in the rriQ —>■ oo. 
Starting from the calculations in the leading order, the decay amplitudes for finite rriQ can 
be calculated unambiguously order by order in the l/mq expansion in HQET. 

Following [^] the local interpolating current for creating 0~ pseudoscalar B meson is 
taken as 

JI-,1/2 = ^ h^l^q , (15) 

and the local interpolating currents for creating 1"*" and 2+ (Di, D2) mesons in the doublet 
(Di, D^) are taken as 



JlTZ = \h drvr + irv^^ - - gT'-'Tpt ) g , (17) 



,3/2 - Y ^ ^-^^ {^t'^T + iT'^T - \ 

where D is the covariant derivative and ^t=l^~ f'^^- Note that, without the last term 
in the bracket in (|T6|) the current would couple also to the 1"'" state in the doublet (0^, 1^) 
even in the limit of infinite mg. 

The QCD sum rule calculations for the correlators of two heavy-light currents give: 

MM 



where ml {qq) = {qgo-^yC^^ q) . 

The QCD sum rule analysis for the sub leading form factors proceeds along the same 
lines as that for the leading order Isgur-Wise function. For the determination of the 
form factor T^^g, which relates to the insertion of Aqcd/'^c kinetic operator of the HQET 
Lagrangian, one studies the analytic properties of the three-point correlators 



I'j d'xd'x'd'ze^^''-'-'-^ {0\T {Jl^,,/,ix') Ol^lA^) J{^^r'\0) 4,-.i/2W) |0) 

= E{uj,uj',y) C^"^^ , (20a) 

^'| d'xdS'd'ze^^''-'-'-^ {0\T {Jf^,,/,ix') O^^IA^) j{^^r'\o) JoV.i/2(^)) |0) 

= S(.;,c.',y)£^"i, (20b) 
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where Jy^^'^ ■* = h{v')'-f'^ h{v), J^^'^ -* = hiy')^^^^ h{v). The variables k, k' denote residual 
"off-shell" momenta which are related to the momenta P of the heavy quark in the initial 
state and P' in the final state hy k = P — rnqv, k' = P' — mgiv' , respectively. For 
heavy quarks in bound states they are typically of order Aqcd and remain finite in the 
heavy quark limit. Cv,a are Lorentz structures associated with the vector and axial vector 
currents(see Appendix). 

The coefficient '^{uj^uo'^y) in (pOD is an analytic function in the "off-shell energies" 
uj = 2v ■ k and u' = 2v' ■ k' with discontinuities for positive values of these variables. 
It furthermore depends on the velocity transfer y = v ■ v', which is fixed in its physical 



region for the process under consideration. By saturating (pO]) with physical intermediate 
states in HQET, one can isolate the contribution of interest as the one having poles at 
UJ = 2A_ 1/2, 1^' = 2A_|_ 3/2. Notice that the insertions of the kinetic operator not only 
renormalize the leading Isgur-Wise function, but also the meson coupling constants and 
the physical masses of the heavy mesons which define the position of the poles. The correct 
hadronic representation of the correlator is 

^ / / X /-,l/2/+,3/2 / / N , 

(2A_^i/2 - - ?e)(2A+_3/2 - cj' - 2e) V 

, -^+,3/2 



(<^+,3/2 + 77a > ~) ^(^) + higher resonance , (21) 

Zi\.4_^3/2 — UJ — J 



where /pj-^ are constants defined in ([TsD, Apj^ = mpj-^ — mg, and Kp^j^ and G^j^ are 
defined by |0],|2l 



(j, P,j,|Ogjj, P,j,) = i^p,,, , (22a) 

m J og,(x) j;-:;^(o)ij, p, j,) = /p,, c^^^r--^ . (22b) 



Furthermore, as the result of equation (|13|), only one state with = or = 2"*" 
contributes to (pi]) , the other resonance with the same quantum number and different 
ji does not contribute. This would not be true for = 1"*" if the last term in (^) is absent. 

Following the standard QCD sum rule procedure the calculations of E{uj,uj',y) are 
straightforward. In doing this, for simplicity, the residual momentum k is chosen to be 
parallel to v such that k^ = {k- v)vn (and similar for k'). Confining us to the leading order 
of perturbation and the operators with dimension D < 5 in OPE, the relevant Feynman 
diagrams are shown in Fig 1. The perturbative part of the spectral density is 

3 1/ / 

Pp,,,{u, u', y) = 28^2 (i^^)(^2_ 1)5/2 ((2^ - 3)'^' + (%' + 12y - 

+ (6y2 + 18?/ - ?,)Cj'^Cj' - {I2y^ + ISy^ + 9)cDcD'2^ 
xQ{uj)Q{Cj')Q{2yCjCj' -Cj^ -Cj'^) . (23) 
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The QCD sum rule is obtained by equating the phenomenological and theoretical ex- 
pressions for H. In doing this the quark-hadron duality needs to be assumed to model 
the contributions of higher resonance part of Eq. (^). Generally speaking, the duality 
is to simulate the resonance contribution by the perturbative part above some threshold 
energies. In the QCD sum rule analysis for B semileptonic decays into ground state D 
mesons, it is argued by Neubert, Block and Shifman in P JT^JT^ that the perturbative 
and the hadronic spectral densities can not be locally dual to each other, the necessary 
way to restore duality is to integrate the spectral densities over the "off-diagonal" variable 
tD_ = (cj — a)')/2, keeping the "diagonal" variable a)_|_ = (a) + cj')/2 fixed. It is in u;+ that 
the quark-hadron duality is assumed for the integrated spectral densities. We shall use the 
same prescription in the case of B semileptonic decays into excited state D mesons. 

The functions in (p3D imply that in terms of cj^ and tD„ the double discontinuities of 
the correlator are confined to the region —\/y^ — 1/ (1 + \/y'^ — 1/ {1 + y) u)+ 

and a)+ > 0. According to our prescription an isosceles triangle with the base (D+ = Cjc is 
retained in the integration domain of the perturbative term in the sum rule. 

In view of the asymmetry of the problem at hand with respect to the initial and final 
states one may attempt to use an asymmetric triangle in the perturbative integral. How- 
ever, in that case the factor (y^ — 1)^/^ in the denominator of (^) is not canceled after the 
integration so that the Isgur-Wise function or it's derivative will be divergent at y = 1. 
Similar situation occurs for the sum rule of the Isgur-Wise function for transition between 



ground states if a different domain is taken in the perturbative integral [|14 



In order to suppress the contributions of higher resonance states a double Borel trans- 
formation in uj and uj' is performed to both sides of the sum rule, which introduces two 
Borel parameters Ti and T2. For simplicity we shall take the two Borel parameters equal: 
Ti = T2 = 2T. In the following section we shall estimate the changes in the sum rules in 
the case of Ti 7^ T2. 

The non-pert urbative power corrections to the correlators are computed from the dia- 
grams involving the quark and gluon condensates in Fig. l(b)-l(k) in the Fock-Schwinger 
gauge x^A^{x) = 0. We find that the only nonvanishing contribution is the gluon con- 
densate. Note, in particular, the vanishing of the mixed quark-gluon condensate {D = 5) 
resulting from the explicit calculation of the diagram shown in Fig. 1(b). After adding the 
non-perturbative part and making the double Borel transformation one obtains the sum 
rule for ?7ke as follows 

,3/2 "T 



in which the sum rule for r(y) has been derived from the study of the three-point correlator 
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in Ref . ill 



as 



r{y) f- 1/2 /+ 3/2 e-(^-.V2+A+,3/.)/T ^ r rf^^ 4 e--+/^ - -m^M 

vy;j ,1/2 J +,3/2 27r2(2/ + l)3 io + 12 ° T 

' ■(-GG)/±4. (25) 



3x25^7r (2/ + 1)2 

From the consideration of symmetry, the sum rule for rjl^'J that originates from the 
insertion of AqcD/^^fe kinetic operator of the HQET Lagrangian is of the same form as 
in (p^), but with the HQET parameters ^^+,3/2 replaced by -^-,1/2, 

respectively. The definitions of and -ft'-, 1/2 can be found in Eq. (|2^). 

It is worth noting that the QCD 0(«s) corrections have not been included in the sum 
rule calculations. However, the Isgur-Wise function obtained from the QCD sum rule 
actually is a ratio of the three-point correlator to the two-point correlator results. While 
both of these correlators are subject to large perturbative QCD corrections, it is expected 
that their ratio is not much affected by these corrections because of cancellation. This has 
been proved to be true in the analysis of Ref. . 

We now turn to the QCD sum rule calculations of the functions parameterizing the time- 
ordered products of the chromomagnetic term in the HQET Lagrangian with the leading 
order currents, rji (z = 1, 2, 3). To obtain QCD sum rules for these universal functions one 
starts from three-point correlators 

^'J d'xd'x'd'ze^^''-'-'-^ (0|T(jr,+,3/2(^') O^^A^) J{^^r'\o) JoV,i/2(^)) |0) 

= -v,A(^^^'^y)^ (26a) 

^'| d'xd'x'd'ze^^''-'-'-^ (o|r(j,"^_3/,(x') o%^AA J^':/\^) ^V,i/2(^)) lo) 

= S^f^^(.;,^',y). (26b) 

By saturating the double dispersion integral for the three-point functions in (p6|) with 
hadronic states, one can isolate the contributions from the double pole at = 2A_. 1/2, oj' = 
2A+ 3/2. Similarly, the insertions of the chromomagnetic operator result in the corrections 
to the leading Isgur-Wise function as well as to the couplings of the heavy mesons to the 
interpolating currents and to the physical meson masses. It follows from Eq. (EBI) that 



'^:fii^polef^ \ _ /-,l/2/ +,3/2 ^ _ ^ 

^ ' ^ ^ (2A_,i/2-^-«e)(2A+,3/2-^'-^e) 

i,C^\ + (G^,3/2 + . ) rfi,3/2r(i/) £r ) , (27) 



^Aii/polo / / N _ /-,l/2/+,3/2 / 

' " (2A_.i/2 - c - ^e) (2A^.3/2 -UJ'- ze) I 



+,3/2 

,3/2 



(2A_,i/2 -UJ- ie)(2A+,3/2 - uj' - 
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— /i"/3pole/ / \ /-,l/2/+,3/2 / ^ 

(2A_,i/2 -UJ- 2e)(2A+,3/2 - u;' - ze) 



+ (G?,3/2 + ^. ) 4,3/2r(l/) ) , (29) 



(2A_^i/2 - u - le) (2A+^3/2 - u' - it) ^ 

-,3/2 



(G^,3/2 + TH^ ■ ) ^2,3/2r(y) , (30) 

^' ' 2A4-.3/2 - cj' - / 



where ^1 = 2?7i + 8773, ^2 = — 16771 — 4(y — 1)772 — 4?73, C, = 2?7i — r73, the quantities Sp^^-^ and 
G'5?p,i, are defined by 

(0|z J d'x O(,«i,(x)j;:^;,7(0)|j,P,j,) = d^/p,,, G^,,,r/°^-"^ , (31b) 
= djj^, dj^^i/2j, = 2ji + 2, dj^+i/2,j, = -2ji , 

Cy", Cy'^, Cy^^ , C^^ , C^^^ dJid C^Ar^ a^^^ defined in Appendix 0. The three-point correla- 
tors ( p6l) can be expressed in QCD in terms of a perturbative part and nonperturbative 
contributions, which are related to the theoretical calculation in HQET. When we do not 
consider radiative corrections, the insertions of the chromomagnetic operator only con- 
tribute to diagrams involving gluon condensates and do not contribute to the perturbative 
diagrams since there is no way to contract the gluon contained in Omag- That is, the 
leading nonperturbative contributions are proportional to the gluon condensates, while the 
leading perturbative contributions are of order and come from the two-loop radiative 
corrections to the quark loop. In general, the calculation of the two-loop diagrams is rather 
cumbersome. In this article we shall neglect the perturbative term of order as and only 
include the nonperturbative gluon condensates without radiative corrections in the sum 
rules. 

Within this approximation one can perform the calculation conveniently by using the 
Fock-Schwinger gauge. After making the double Borel transformation, the sum rules for 
rji{y) are obtained as follows 

r?2(y)=0, (32b) 

3/2 "I ^jT 

^ ■(^GG)-ii^T . (32c) 



/-,i/2 /+,3/2 e-(^-.v^+A,,3/.)/T ^ 2(G^,3/2 + ^)r(y) /_,V2 /+,3/2e-(^->v^+A,,3/.)/T 



240 ' TT ' {y + 1) 
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^71,2,3 are expected to be small compared to Aqcd since the mass splitting between and 
Di is very small. This is supported by the fact that the QCD sum rule calculations indi- 
cate that the analogous functions parameterizing the contributons of the chromomagnetic 
operator for B D^*^ eue decays are small [0. The Aqcd /nib correction associated with 
the insertion of chromomagnetic operator of the HQET Lagrangian can be investigated in 
a similar way. 



IV. NUMERICAL RESULTS AND IMPLICATIONS FOR B DECAYS 



We now turn to the numerical evaluation of these sum rules and the phenomenological 
implications. For the QCD parameters entering the theoretical expressions, we take the 
standard values P,PD| 



{qq) = -(0.23 ± 0.02)^ GeV^ , 

{^GG) = (0.012 ± 0.004) GeV^ , 
vr 

ml = (0.8 ± 0.2) GeV^ . (33) 

In order to obtain numerical results for r^ked/), Vkdv) rii{y) {i = 1,2,3) from the sum 
rules which are independent of specific input values of /'s, A's and r, we adopt the strategy 
to evaluate the sum rules by eliminating the explicit dependence on these quantities by 
using the sum rules for them. Substituting the sum rules ([T8|) and (|1^) into the left side 
and the sum rule (|25|) into the right side of the sum rules (|2^) and (|3^) for the three-point 
correlators, we obtain expressions for the T^ke, Vke rji {i = 1,2,3) as functions of the 
Borel parameter T and the continuum thresholds. This procedure may help to reduce the 
uncertainties in the calculation. For other HQET parameters we use the following values 
obtained by QCD sum rules 



ir+,3/2 = -(2.0 ± 0.4) GeV^ , G^^y^ = -(1-0 ± 0.45) GeV 

Js:_,i/2 = -(1.2 ± 0.20) GeV^ , = -(1-6 ± 0.6) GeV , 

= (0.020 ± 0.003) GeV^ , G+^g/a = (0.013 ± 0.007) GeV , 

E_,i/2 = (0.23 ± 0.07) GeV^ , G^ ^^ = (0-042 ± 0.034 ± 0.053) GeV . (34) 

Let us evaluates numerically the sum rule for rjy^eiy) and rj^civ) at first. The continuum 
thresholds Ud and Uc2 in ( |T8l) and (|19D are determined by requiring stability of these sum 



rules. One finds that 1.7 GeV < lu^i < 2.2 GeV and 2.7 GeV < < 3.2 GeV |T|,|0 
Imposing usual criterion on the ratio of contribution of the higher-order power corrections 
and that of the continuum, we find that for the central values of the condensates and 
HQET parameters given in ( p3D and (|34D, if the threshold parameter Uc lies in the range 
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1.9 < LUc < 2.5 GeV, there is an acceptable "stability window" T = 0.8 — 1.0 GeV in which 
the calculation results do not change appreciably. This window overlaps largely with 
those of the sum rules (|18|) and (|T^) . Therefore, our procedure of calculation is justified. 
For estimating the errors induced by the uncertainties of parameters for the condensates 
and HQET (|33|) and (|3^) we take the maxima deviations from the central values of the 



condensates and HQET parameters and find that for the existence of stability windows 
in the two extreme cases the continuum thresholds shift to the range 2.5 < ujc < 2.9 and 
1.5 < uJc < 1.9 GeV, respectively. The corresponding windows for Borel parameters are 
1.0 < T < 1.2 and 0.6 < T < 0.8 GeV, respectively. These are still compatible with the 
stability windows for the sum rules ([T8|) and (|1^). 

The numerical results of the form factors f]kc{y) and ri^f,{y) are shown in Fig. 2, where 
the curves refer to various choices for the continuum thresholds and to the central values 
of the condensates and HQET parameters. 

The numerical analysis shows that rjkciy) is a slowly varying function in the allowed 
kinematic range for B — > Diiu and B —>■ D^^^v decays. The resulting curve for rjy^eiy) inay 
be well parameterized by the linear approximation 

vUy) = %e(l) (1 - pjd/ - 1)) , %e(l) = 0.38 ± 0.17 GeV , = 0.8 ± 0.1 , 

vLiy) = ^ke(i) (1 - p'Ay - 1)) , ^ke(i) = o.48 ± 0.21 GeV , p^, = 1.0 ± 0.1 . (35) 



The final sum rules for rii(y) can be obtained by substituting Eq. (^) into Eq. (32a) 
and ( |32c| ). The numerical evaluation for these sum rules proceeds along the same lines 
as that for 77kc(y)- Note that we have not included the perturbative term of order as, 
which is the leading perturbative contribution for rii{y). The sum rules for rii{y) are not 
quantitatively reliable. Nevertheless, they are of correct order of magnitude. The values 
of the form factors rii{y) and rj^ly) at zero recoil as functions of the Borel parameter are 
shown in Fig. 3, for three different values of the continuum threshold Uc- The numerical 
results for rii{y) and rj^ly) at zero recoil in the working regions read 

r^i(l) = -0.95 X 10"% 773(1) = 3.5 X 10~2 (36) 

This result is in agreement with the expectation based on HQET that the spin-symmetry 
violating corrections described by T]i{y) are negligibly small. 

Using the forms of linear approximations for f]ke{y) together with T{y) and Ti^2{y) given 
in Ref. iTBS 



T{y) = 0.74(1 - 0.9{y - 1)) , n{y) = -0.4(1 - lA{y - 1)) , T2 = 0.28(1 - 0.5(y - 1)) (37) 

and neglecting the contribution of chromomagnetic correction, we can calculate the total 
semileptonic rates and decay branching ratios by integrating Eqs. (0) and (^. We use 
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the physical masses, ms = 5.279, raD^ = 2.422 and raj:,* = 2.459 ||3^, for B, Di and 
mesons. The maximal values of y in the present case are y^ax = (1 + ^i)/2'"i ~ 1-32 and 
ymax = (1 + r|)/2r2 ~ 1.31. The quark masses are taken to be = 4.8 GeV, nic = 1.5 
GeV. In Table | we present our results for decay rates both in the infinitely heavy quark 
limit and taking account of the first order 1 / mq corrections as well as their ratio 

From Table | we see that the B Diiu decay rate receives large l/mg contributions 
and gets a sharp increase, while the B — > Dg^i/ decay rate is only moderately increased by 
subleading l/mg corrections. The reason for this is as following. From Eqs. (j^), and 
(|l^) we see that {y — ri)'^fy^ term dominates the differential width for decay to Di near 
y = 1. fvi vanishes at the leading order and receives non- vanishing contributions from first 
order heavy quark mass corrections: 

v^/y,(l) = -8£,(A'-A)r(l). (39) 

Since the allowed kinematic ranges for B — > Diiu is fairly small, the contribution to the 
decay rate of the l/mg corrections is substantially increased. On the other hand, the 



matrix elements ([TcD and ([THD of the B D^^dv decay vanish at zero recoil without using 
the heavy mass limit. The term [y — r2Yk\^ dominates the B — > 02^^ decay rate, but /c^i 
does not vanish at the leading order recoil. Therefore, this process is not much affected 
by next-to-leading corrections. Note that although the correction to the rate for decay to 
Di is very large it comes mainly from the effect of the different masses of B and Di. The 
values of r]ke and T^^g in (23) are of the order Aqcd and perfectly normal. 

In Table | the available experimental data for semileptonic B decay to excited D** 
mesons are presented. As for the B Dg^z/ branching ratio there are only upper limits 
from these experimental groups except the data from OPAL. In comparison with the exper- 
imental data our result for the branching ratio of the B Diiu decay with the inclusion 
of l/mq corrections is larger than the CLEO and ALEPH measurements but is consistent 
with OPAL and DELPHI data. On the other hand, our branching ratio for the B D2ii^ 
decay disagrees with the ALEPH data but is consistent with results from other groups. 



V. CONCLUSION 

In this work we have presented the investigation for semileptonic B decays into excited 
charmed mesons. Within the framework of HQET we have applied the QCD sum rules to 
calculate the universal Isgur-Wise functions up to the subleading order of the heavy quark 
expansion. The Isgur-Wise functions ?7ke and rj^^ related to the insertions of kinetic energy 
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operators of the HQET Lagrangian are found of normal values of the order Aqcd, while the 
form factors rji, parameterizing the time-ordered products of the chromomagnetic operator 
in the HQET Lagrangian with the leading order currents, are neghgibly small.These results 
are in agreement with the HQET-based expectations. 

We have computed, for the decays B —>■ Diiu and B — > D^^f-u, the differential decay 
widths and the branching ratios with the inclusion of the order of l/mq corrections. Our 
numerical results show that the first order l/mq correction is not large for the decay rate 
oi B D2iP process, but is very large for the B DiiP process. We have explained the 
reason for this result. 
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APPENDIX A: 



We list here the lorentz structures used in the paper. 

1 



l-'V 



{y^ - l)gr + (3^;'^ + (2 - y)v'^y, 



+ ^"^ {v'ivt-\{\-y')gf 



C^Tl = (-^^ + v'n (2<^f - ^(1 - y')9f) , 

where gf = g"!^ - v"^v'^ and ^ v" - yv"^. 



(Al) 
(A2) 
(A3) 

(A4) 
(A5) 
(A6) 
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Figure Captions 



Fig. 1. Feynman diagrams contributing to the sum rules for the Isgur- 
Wise form factor in the coordinate gauge. The gray square corresponds 
to the insertion of the kinetic energy operator at 0{l/mQ) in the HQET 
Lagrangian. 

Fig. 2. Results of the numerical evaluation for the sum rules: Isgur- 
Wise form factors r?ke(2/) ^^id r^^{y) with T = 0.9 GeV. 

Fig. 3. Dependence of r]i{l) and 772(1) on the Borel parameter T for 
different values of the continuum threshold u>c- 
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TABLES 

TABLE I. Decay rates T (in 10"^^ GeV) for \Vcb\ = 0.04 and branching 
ratios BR (in % and taking tb = 1.6ps) for B D**£u decays in the infinitely 
heavy quark mass hmit and taking account of first order l/mg corrections. 







B Diev 


B D^ev 


ruQ oo 


r 


1.4 


2.1 




Br 


0.34 


0.52 


With l/niQ 


r 


5.3 


2.4 




Br 


1.3 


0.59 






3.78 


1.15 


Experiment 


Br (CLEO) 


0.56 ±0.13 ±0.08 ±0.04 


< 0.8 




Br (ALEPH) § 


0.74 ±0.16 


< 0.2 




Br (OPAL) § 


2.0 ±0.6 ±0.5 


0.88 ±0.35 ±0.17 




Br (DELPHI) | 


] 1.5 ±0.55 


< 6.25 
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